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Abstract. 

Measuring correlations between forward dijets in deuteron-gold collisions at RHIC 

will further test the Color Glass Condensate (CGC) picture of the nuclear wavefunction 

at small— x, which successfully predicted the suppressed production of high-p<r particles 

■ \ at forward rapidities in d+Au collisions: BdA < 1- I present predictions for the 

*j ■ correlation in azimuthal angle between forward dijets in the CGC framework, with 

^S^ , both multiple scatterings and non-linear QCD evolution at small— x taken into account. 

>, 

£1? ■ 1. Forward particle production in p+A type collisions 

00 ' 

Forward particle production in proton-nucleus collisions allows one to investigate the 

non linear QCD dynamics of high-energy nuclei with a probe well understood in QCD. 

Indeed, while such processes are probing small-momentum (small— x) partons in the 

nuclear wavefunction, only high-momentum partons of the proton contribute to the 

scattering. For instance, in the case of two-particle production, the longitudinal fraction 

of momenta probed in the proton (x p ) and in the nucleus (xa) are given by 

^x p = \ki\ e yi + \k 2 \ e y2 , v^x A = N e- yi + \k 2 \ e~ m (1) 

where kx, k 2 and y\, y 2 are the transverse momenta and rapidities of the final state 
particles. Therefore, with ^/s ^> \ki\, \k 2 \ and forward rapities y%,y 2 > 0, the process 
features x p < 1 and i^Cl, meaning that the scattering involves a well-known dilute 
hadron p and dense nuclear target A, whose non-linear QCD dynamics can be studied. 

The Color Glass Condensate (CGC) framework was quite successful in describing 
single inclusive particle production at forward rapidities in d+Au collisions at RHIC 
[H [2]. In this work, we focus on forward inclusive two-particle production pA^hih 2 X 
with h\ and h 2 detected in the proton direction, and in particular on correlations in 
azimuthal angle between the produced hadrons |3j. We argue that the second d+Au 
run at RHIC gives the opportunity to carry out these measurements, and check the 
relevance of the CGC picture of high-energy nuclei at RHIC energies. 

In the following, the pA^hih 2 X cross-section is derived and used in the context 
of d+Au collisions to predict the azimuthal angle distribution. 
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2. Forward dijet production 

The kinematic ranges for forward particle measurements at RHIC are such that a; p ~0.5 
and xa ~ 10~ 4 . Therefore the dominant partonic subprocess is initiated by valence quarks 
in the proton and, at lowest order in a s , the pA^hiti2X cross-section is obtained from 
the qA — > qgX cross-section, the quark density in the proton and the appropriate hadron 
fragmentation functions. 

By contrast, the nucleus A, whose partons with small fraction of momentum (mainly 
gluons) are described by a CGC, cannot be described by a single gluon density The 
qA — > qgX cross section is instead expressed in terms of correlators of Wilson lines 
(which account for multiple scatterings), with up to a six-point correlator [4j. Assuming 
a Gaussian CGC wavefunction allows to express all the correlators in terms of a single 
function T(r, xa), related to the variance of the Gaussian distribution [5]. In the 
large— N c limit, the non-linear QCD evolution at small— xa is included through the 
Balitsky-Kovchegov (BK) evolution P,[7j of e~ r . 

The cross-section for the production of the quark-gluon dijet (with respective 
transverse momenta q± and k± and rapidities y q >0 and yk>0) then reads [4]: 

J^pA—tqgX ry C 1 

JM^dwto, * -^ x ' q{x " A) g ^o {z ' kj -~ A; **)-*&'■ k -~^)\ F »( A > ( 2 > 

with A = kj_ + gj_ and 1/z — 1 = \q±\ e Vq ~ Vk /\k±\. The different components in ([2]) are: 

• the quark distribution function in the proton q(x p , A) (we are working with 
|A| 3> Aqcd)'i i n principle the gluon initiated processes gA — ► qqX and gA — > ggX 
should also be included, they would contribute for smaller values of x p . 

• the Fourier transform of the two-point (dipole) correlator F XA (A), also called the 
unintegrated gluon distribution function, given by 

F - (A) = /(2^ e " iA ' re_r(r '" ); (3) 

contrary to the dilute proton, the CGC cannot be described only by its gluon 
distribution F XA (A), the kx— factorization formula is not applicable 0, [9] (it would 
if there was no other xa dependence in (J2J), as in single gluon production [TC)l[Tl"] ). 

• the kx~ factorization breaking terms (whose simple form is due the use of a Gaussian 
CGC wavefunction) \. \^ap{ x A) — ipap] w ith 

Act/3 

I^(z,k±;x A ) = / d 2 q ± ip* f3 (z,q ± )F XA (k±-q±) ; (4) 

as the valence quark emits the virtual gluon (the associated wavefunction is ip^ 
where A, a and f3 are polarization and spin indices) it interacts coherently with the 
dense small— x gluons in A, which modifies the q^qg splitting. 

In practice, the initial condition F(r, Xq) = r 2 Q 2 So ln[e + l/(r 2 AQ CD )]/4 is used with the 
choice Xq = 0.01 and the initial saturation scale 2-n Q 2 SQ = 2 GeV 2 . Then e _r is evolved 
with the BK equation, allowing to compute (J3J), (jlj) and the cross-section fl2]). 
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3. Dijet correlations in azimuthal angle 

We will now use the inclusive two-particle spectrum (J2J) to investigate the process 
dAu — > hih,2X, with yfs = 200 GeV. In particular we shall study the A0 spectrum 
where A0 = 0i — 02 is the difference between the azimuthal angles of the measured 
particles hi and h?.. We will study the normalized A0 distribution 

1 da / da hT^ hlh2 x \ -i d(J hT^h,h 2 X 

a dAcj) \dpT 1 dpT 2 dyidi/2 J dpT 1 dpT 2 dyidy2dA(j) 

where k\ = (pti,0i) an d &2 = (PT 2 ,(p2) are the transverse momenta of the measured 
hadrons and ?/i and y2 are their rapidities. We take into account both situations where 
either the quark or the gluon fragments into the most forward hadron. Our results do 
not include convolution with fragmentation functions which should be done eventually, 
however we except this has little impact on the observable we are considering here; for 
Pt~ spectra this would be not be case. 

The main features of the A0 spectrum (JSj) are the following. First we obtain that 
the perturbative back-to-back peak of the azimuthal angle distribution is recovered for 
very large transverse momenta. Then, this back-to-back correlation is reduced by the 
initial state saturation effects included in our CGC approach: as the transverse momenta 
decrease closer to the saturation scale (Q s — 2 GeV), the angular distribution broadens. 
Finally, we notice that at RHIC energies, this does not lead to a complete disappearance 
of the back-to-back peak. 

As can be seen from the kinematics ([IT), the most forward of the two particles 
essentially determines the value of Xh while the most central one determines the value of 
xa- In order to quantitatively study the effect of the CGC evolution, the ideal situation 
would be to keep Xh fixed and to vary xa- In practice, this is better realized by fixing the 
rapidity and momentum of the most forward particle and by varying the kinematics of 
the other. Note that doing the opposite would emphasize the Xh evolution of q(xh, A), 
rather than focus on the xa evolution of F XA . Moreover, the cross-section (J2J) is quite 
sensitive to choice of factorization scale in the quark density, so it is better to keep 
Xh constant. Note that varying the rapidities at fixed y\ — y% would keep the product 
XhXA constant, and would force a competition between the evolution of q(xh, A) with 
increasing Xh and the CGC evolution with decreasing xa- 

In Figure la, we have studied the A0 spectrum (jSJ) in the situation in which 
Pt-! = 3.5 GeV, px 2 = 2 GeV, y\ = 3.5 and y 2 is varied from 1.5 to 2.5. As y 2 increases, 
the value of xa decreases and the suppression of the azimuthal correlation is more 
important. However the effect is quite small, because the increase of the saturation 
scale with decreasing xa is rather slow. In Figure lb, we investigate the situation for 
which px 1 =5 GeV, y\ = 3.5, y 2 = 2 and p? 2 is varied 1.5 GeV to 3 GeV. As pt 2 decreases, 
it gets closer to the saturation scale Q s (which also slightly increases as xa decreases), 
and the suppression of the azimuthal correlation increases. Varying pt 2 at fixed y^ allows 
to probe the ratio Pt 2 /Qs over a larger range, so the effect is much bigger than when 
varying y 2 at fixed pr 2 ■ 
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pr,=3.5 GcV, pt,=2 GcV, yi = 3.5 




PT,=5GeV, yi = 3.5, y2-2 



pr 2 =3.0GeV - 
pii=2.5 GeV - 
PT a = 2.0GeV - 
pr„ = 1.5 GeV ■ 




Figure 1. The Acf> spectrum ([5]) in two situations with the RHIC energy y/s = 
200 GeV/nucleon. In both cases, by increasing yi or decreasing pt 2 , the correlation 
in azimuthal angle is suppressed as the value of xa probed in the process decreases. 
Varying pt 2 at fixed 2/2 is much more efficient as the ratio pr 2 /Qs covers a larger range. 



Experimental measurements of two-particle correlations in azimuthal angle have 
been performed in d+Au collisions at RHIC by the PHENIX and STAR collaborations. 
The measurements of STAR [12] with tt° at forward rapidity and charged hadrons at 
mid rapidity are qualitatively consistent with a suppression of the back-to-back peak 
with respect to p+p collisions. By contrast, the measurements of PHENIX [13] do not 
show any evidence of a suppression of the back-to-back peak, but they probe values of 
xa which are bigger than 0.01. It may very well be that the CGC picture breaks down 
for values of xa bigger than 0.01, and it justifies our choice not to start the small— xa 
evolution at a higher value. 

Future measurements of azimuthal correlations between forward particles in d+Au 
collisions, which could be carried out at RHIC, will allow quantitative tests of the CGC. 
Our predictions for the fully differential cross section ([2]) are not directly comparable 
with data: once the cuts used by the experiments become available, integrations over 
the kinematic variables should be performed. Measurements in p+Pb collisions at the 
LHC would reach x^^lO -5 , and could test even better the QCD evolution at small— x. 
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